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Abstract 

We describe in detail a general strategy for implementing a conditional 
geometric phase between two spins. Combined with single-spin opera- 
tions, this simple operation is a universal gate for quantum computation, 
in that any unitary transformation can be implemented with arbitrary 
precision using only single-spin operations and conditional phase shifts. 
Thus quantum geometrical phases can form the basis of any quantum 
computation. Moreover, as the induced conditional phase depends only 
on the geometry of the paths executed by the spins it is resilient to certain 
types of errors and offers the potential of a naturally fault-tolerant way 
of performing quantum computation. 

1 Introduction 

Among the surprising effects recently discovered in quantum mechanics is that 
a quantum system retains a memory of its motion when it undergoes a cyclic 
evolution This is reflected in the existence of the Berry phase, a phase 
acquired by the quantum state of the system in addition to the better known 
dynamic phase. The Berry phase is a purely geometrical effect that can be linked 
to the notion of parallel transport it depends only on the area covered by 
the motion of the system, and is independent of details of how the motion is 
executed. Berry phases have been demonstrated in a wide variety of systems 
[D, including NMR ||,||, the closely related technique of NQR optical 
systems and others. 
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An equally exciting recent development in the field of quantum mechanics 
has been the discovery that quantum systems can be used to perform novel 
information processing tasks, including computations which are more efficient 
than any algorithm known on a classical computer pO[-p^. Quantum informa- 
tion processing requires the ability to execute conditional dynamics between 
two quantum bits (qubits), where the state of one qubit influences the evolution 
of another qubit during a quantum computation. Simple quantum information 
processing has been demonstrated using NMR [p^-p7| and trapped ions |p^ . 

Recent experimental work has managed to combine these two quantum phe- 
nomena in the form of geometric quantum computation ||l9[| . (For a more ab- 
stract approach see |^,^.) In this paper we seek to detail the theoretical ideas 
behind geometric quantum computation. In particular we show that Berry's 
phase may be used to implement conditional phase shifts, and thus any quan- 
tum gate HJ] . We begin with brief introductions to both quantum gates and 
networks as well as to geometric phases, proceeding to analyse the dynamics of 
a spin-half system in order to see in detail how the theory of geometric phases 
applies there. Finally, we extend the ideas to pairs of spin- half particles, showing 
how to introduce a conditional geometric phase between the two particles. 



2 Phase gates and quantum computation 

2.1 Qubits and networks 

A qubit is a quantum system in which the Boolean states and 1 are represented 
by a prescribed pair of normalised and mutually orthogonal quantum states 
labeled as {|0),|1)}. Unlike a simple Boolean variable, a qubit, typically a 
microscopic system such as an atom, a nuclear spin, or a polarised photon, can 
exist in an arbitrary superposition a | 0) + /3 | 1), making it more powerful as a 
computational resource. 

In quantum computation, we set some register of qubits to an "input" state, 
evolve the qubits unitarily using simple building-block operations and then take 
the final state as "output" . More formally, a quantum logic gate is a device 
which performs a fixed unitary operation on selected qubits in a fixed period of 
time and a quantum network is a device consisting of quantum logic gates whose 
computational steps are synchronised in time | p3[ |. The outputs of some of the 
gates are connected by wires to the inputs of others. The size of the network is 
the number of gates it contains. 



2.2 Quantum logic gates 

The most common quantum gate is the Hadamard gate, a single qubit gate H 
performing the unitary transformation known as the Hadamard transform. It 
is defined as 
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The matrix is written in the computational basis {| 0) , | 1)} and the diagram on 
the right provides a schematic representation of the gate H acting on a qubit 
in state | x), with a; = 0, 1. 

The addition of another single qubit gate, the phase shift gate 0, defined as 
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1), or, in matrix notation, 
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is actually sufficient to construct the following network (of size four), which 
generates the most general pure state of a single qubit (up to a global phase). 
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Consequently, the Hadamard and phase gates are sufficient to construct any 
unitary operation on a single qubit. 

Thus the Hadamard gates and the phase gates can be used to transform the 
input state |0)|0)...|0) of n qubits into any state of the type l^*!) |^'2)--- |*n)) 
where is an arbitrary superposition of |0) and |1). These are rather special 
n-qubit states, called the product states or the separable states. In general, a 
register of n qubits can be prepared in states which are not separable, known 
as entangled states. 

However, in order to entangle two or more qubits it is necessary to have 
access to two-qubit gates. One such gate is the controlled phase shift gate i?((/)) 
defined as 
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The matrix is written in the computational basis {| 00) , | 01) , 1 10) , 1 11)} and 
the diagram on the right shows the structure of the gate. 



2.3 Universality 

An important result in the theory of quantum computation states that the 
Hadamard gate, and all B{(j)) controlled phase gates form a universal set of 
gates: if the Hadamard gate as well as all B{(l)) gates are available then any 
n-qubit unitary operation can be simulated exactly with less than C4"n such 
gates, for some constant C |Q. Consequently, being able to implement 1 and 
2-qubit phase gates is of crucial importance in quantum computation. In this 
paper we describe a new method for implementing the controlled phase gates 
based explicitly on geometric phases rather than dynamic ones. 
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3 Geometric phase 



3.1 Cyclic evolution 

The states of a quantum system are usually described as being represented by 
vectors of norm 1 (| ( ■0 | ) p = 1) in a complex Hilbert space Ti. However, there 
is redundancy in this description since the state | ifi) is physically indistinguish- 
able from the state e^'^ | "i/i) . It is therefore convenient to consider instead the 
projective space 7-", in which vectors are grouped into equivalence classes under 
the relation | ■0) ~ re**^ | ip) for any r > and real (/>, thereby eliminating the 
ambiguity. The associated projection map is 



n 



n 

IV') 
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(5) 



If a system undergoes a cyclic evolution, the ket representing the system 
state traces out a path, C : [0, r] — > 7i, where 11 (C) is a closed curve in T', as 
illustrated in Figure [|. In other words, the initial and final states should be on 
the same ray in 7i, but may be related by a phase, e**^. We will measure this 
phase with respect to a reference curve in Ti.: for each point | on C, we can 



choose a smoothly varying representative 
that -0(0)) = tl!{T)). We can then write 



il^{t)) from Il{'tp{t)) in such a way 



(6) 



so that the phase change of | ■0(0)) associated with the cyclic evolution is given 
by0 = /(r)-/(O). 

3.2 Dynamic and geometric phase 

The time evolution of a quantum system is governed by the Schodinger equation. 



in-m)) = H{t)\m). 



(7) 



where H{t) is the Hamiltonian. Substituting Eq. (^ into the above, rearranging 
and multiplying by {ip{t) \ gives the following 123], 



'^ = -\{m\H\m)+^{m 
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or, when integrated 
1 



^^j^ m)\H\m)dt+ij^ ^ mjdt. (9) 



Thus, (j) can be decomposed into a dynamical phase 

^ ''\m\H\m)dt 



(10) 
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Figure 1: A schematic diagram of a spin-half particle undergoing a cyclic state 
evolution. Points on the sphere correspond to physically distinguishable states. 
Going through each point is a ray, on which phase information and normalisation 
is recorded. 
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which depends on the Hamiltonian, and a geometric phase 

-f^illjj d (11) 



which depends only on the path C; 7 is independent of the rate at which | ip{t)) 
progresses along C, the Hamiltonian, or the choice of reference | ^)}- 

3.3 Berry's phase 

A particular instance of this geometric phase is Berry's phase [Q , which occurs 
when the adiabatic theorem (see [^) is satisfied. In this case, if the initial state 
I "0(0)) of the system is an eigenstate of the Hamiltonian, the state | ip{t)) remains 
an eigenstate | ip{t)) = \ n(R)) of the instantaneous Hamiltonian -ff (R), where 
R is a set of time- varying parameters controlling the Hamiltonian. Supposing R 
traces a closed loop in the parameter space, the geometric phase of the system 
can be written in terms of R. In this case, provided the energy eigenspace of the 
instantaneous Hamiltonians is non-degenerate along the path C, the geometric 
phase acquired by the n"^-eigenstate is 



-fn^if [n{R) Vr n(R)^ ■ dR, (12) 

where Vr is the gradient operator with respect to the parameters R, and 
n(R)^ is defined as in the previous section. This line integral can be trans- 
formed into a surface integral over any surface in the parameter space whose 
boundary is C. 

Since experimentally it is much easier to control the Hamiltonian than the 
actual state of a system, the adiabatic case is of importance. However, the 
adiabatic conditions necessarily mean that the processes take a long time com- 
pared to the characteristic dynamical time-scales, and thus are much slower 
than dynamic methods of generating phases. 



4 Single-qubit evolution 

4.1 Qubit Dynamics 

Here we will focus on developing an understanding of the time evolution of 
a single qubit governed by a very general Hamiltonian. Recall that any 2x2 
Hermitian matrix can be written in terms of the unit matrix and the three Pauli 
matrices. In particular a single qubit density operator can be parametrised as 

where the real vector s = {sx, Sy, Sz) is known as the Bloch vector. By the same 
token any 2x2 Hamiltonian can be written as 

H='^{nol + n-a), (14) 
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Figure 2: Solution to the equations of motion for a single spin-half particle. 



where ft is called the Rabi vector. Substituting Eq. ( |13| ) and Eq. (|14| ) into the 
equation of motion for the density operator, 

d 



and using the identity 

(a • cr)(b • ct) = (a • b) 1 + i{a xh) ■ a, 
we find the following equation of motion for the Bloch vector, 

d 



dt 



r2 X s. 



(15) 
(16) 
(17) 



This equation has a simple geometric solution: vector s revolves around vector 
r2 with angular frequency given by the length of $7, as illustrated in Figure 

From Eq. (|17D, it is relatively easy to move to situations typical of those 
encountered in quantum computation. Hamiltonians that describe qubits in- 
teracting with external potentials are usually time dependent. Typical exter- 
nal perturbations are periodic such as, for example, spins coupled to oscillating 
magnetic fields in NMR or atomic dipole moments coupled to oscillating electro- 
magnetic field in the optical domain. Within the Rotating Wave Approximation 
the oscillating field can be replaced by a rotating field, and so the Hamiltonian 
is of the form 



Hit) = 



h 



LOie 



(18) 



where ujq /2Tr is the system's transition frequency, while lo /2t: and fkui are the 
frequency and the amplitude of the oscillating field, respectively. This gives 

fix = ij-'iCOs{ujt + (jj), fly = ui sm{ujt + (j)) , flz = loq. (19) 
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In order to solve Eq. ( J17| ) it is convenient to consider the evolution of s in 
a frame which rotates with frequency lu around the z-axis. More precisely, we 
write 

s{t) = R,{ujt)s'{t), n{t) = R,{ujt)n'{t), (20) 
where Rziujt) is the rotation matrix 

(cos(ti)i) — sin(ti)i) \ 

sin(wt) cos(tjt) = exp(tjtM^), (21) 

1 y 

for 

/O -1 0\ 

= 1 . (22) 

\ / 

Substituting Eq. ( pO| ) into Eq. and taking into account that 

i?,(wO = R,{ujt){LoM,), M,s' ^ z X s', (23) 

where z is a unit vector in the z-direction, we obtain 

d 



d^^ 



n' X s' (24) 



with the time-independent vector Jl', 

= wi cos(0), 17^ = sin(0), fi'^ = — uj. (25) 

If we can control the strength of coupling fkoi, the frequency u and the phase 
(j) of the external field we can prepare any vector Jl. This implies that if we 
know the initial state of the qubit then with a single rotation we can position 
the Bloch vector s in any prescribed direction. 

4.2 Calculating geometric phases 

We can now apply the results of the previous section to the task of developing 
a geometric phase of a spin-half particle located in an external oscillating field. 
By varying the parameters of the Hamiltonian adiabatically we will send a qubit 
through a cyclic evolution, whose associated geometric phase can be calculated 
using Eq. (|ll|) . 

Working in the rotating frame, the components (17^,17^,51^) of the Hamil- 



tonian are given by Eq. (25), where the frequency and power of the rotating 
field can be used to control the value of the angle 9 of Figure ||. In the absence 
of a rotating field the Rabi vector lies along the z-axis, and as the power of 
the rotating field is slowly increased from zero to cui the vector tilts towards the 
xy-plane. If the Bloch vector s' is initially aligned with ft' then by the adiabatic 
theorem, it will remain aligned with $7' provided fJ' varies slowly. Therefore, 
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Figure 3: Spin-half particle in a magnetic field 



from the relation between the different components ft' we find that the angle 9 
between the Bloch vector and the ^-axis will be 

cos^^ . ^ . ... (26) 



Varying the phase </> of the rotating field from Eq. (25) will then cause 
s' to rotate around the z-axis. The geometric phase associated with a full 
revolution is easily calculated by parameterizing the state in terms of the 1 1) 

and 1 1) eigenstates of quantization along the z-axis, ?/'(a)^ = cos(0/2) 1 1) + 

sin(0/2)e'" | J,), where a changes smoothly from to 2tt. Using Eq. (ITTl), 



^ £ (^cos ^ I T) + sin ^e^ \ 1)^ d (^cos ^ 1 1) + sin | j) 



cos - (t I + sin -e-» (i | j (^sin -e^ | i) j da 

= --(1 -cos6l) / da 
2 Jo 

= -7r(l-cos6l). (27) 

Furthermore, the result can be generalized to any closed path with the result 
that the geometric phase is equal to half the solid angle enclosed by C on the 
Bloch sphere [Q . 
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4.3 Eliminating dynamic phases 



In order to perform conditional quantum gate operations using geometric phases 
only, it is necessary to find a way to eliminate the dynamic phase. One approach 
is to use a refocussing technique known as spin-echo. The basic idea is to apply 
the cyclic evolution twice, with the second application surrounded by a pair 
of fast TT transformations (this being simply the transformation that swaps the 
basis states 1 1) and 1 The net effect of this compound transformation would 
be to cancel all the acquired phases except that the second cyclic evolution is 
performed by retracing the first but in the opposite direction so that while the 
dynamic phases cancel, the geometric phases do not. 



I ) during 



To see why this is so, let C| be the closed curve in V traced out by 

the first cyclic evolution and the one traced out by J^, with corresponding 

dynamic and geometric phases Si, 7^, Si, and 71 . Referring back to Eq. ( p7| ) 
we see that 71 = —7 and 7x = 7 for 7 = 7r(l — cos 9). Similarly, if we write C-\ 
and Ci for the second cyclic evolution, these are simply C| and C| carried out in 
opposite orientations so they have corresponding phases 7f = 7 and Ji ~ —7. 

In summary, we can follow the states through the compound evolution as 
follows: 

I e^iST--,) I e»(<5T-7) I I) A> gj(<5T +5^-27) I 1^ ^i{S-,+Si-2t) | 

I i)-^ e*('^i+'^) I i)-^ e''(*i+^) I t)-^ g»(<5T+'5i+27) I gi(5T+<5i + 27) | |^ _ 

(28) 

Since the global phase factor e'^'^T+'^i) is not physical this sequence of operations 
behaves as promised. The dynamic phases are eliminated and we are left with 
an exclusively geometric phase difference of 47 = 47rcos6'. 



5 Conditional dynamics 
5.1 2-Spin Hamiltonian 

This geometric phase can be used to implement a 2-qubit controlled-phase gate. 
Consider to begin with a system of two non-interacting spin-half particles Sa 
and Sb. In a reference frame aligned with the static field, the Hamiltonian reads 



Hq = hUaSaz Ife + fK^b^a <^ Sbz, 

or, in the basis {| Sa., Sbz)}s^^,s,. - {I TT) , I Ti) , lit) , I ii)}, 

/ UJa+i^b 

— Wfc 

-UJa +UJb 



V 











(29) 



(30) 



/ 



where the frequencies LOa/lir and ujb/lTT are the transition frequencies of the two 
spins and we have used the scaled Pauli operators Si = (Ti/2. (From now on we 
assume that uJa and LOb are very different with Ua > c^b-) 
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Figure 4: The energy diagram of two interacting spin-half nuclei. The transition 
frequency of the first spin depends on the state of the second spin. 



If the two particles are sufficiently close to each other, they will interact, 
creating additional splittings between the energy levels. In the case of two spin- 
half particles, the magnetic field of one spin may directly or indirectly affect the 
energy levels of the other spin; the energy of the system is increased by TrhJ/2 if 
the spins are parallel and decreased by nhJ/2 if the spins are antiparallel. The 
Hamiltonian of the system taking into account this interaction reads 



H ^ Ho + 2-KhJSaz ® Sbz, 

or, in the previously chosen basis. 



(31) 



H 



\ 





LUa — LOh — TtJ 

—UJa + ^6 — T^J 

— CJa ~ i^b + T^J J 



(32) 



Figure ^ illustrates the energy levels of the system. When spin St is in state 
1 1), the transition frequency of the spin Sa is 

LU^ = LUa + J, (33) 

whereas when spin Si, is in state | J,), the transition frequency of the spin Sa is 

LU^ ~ LUa — TT J. (34) 

5.2 Conditional phase shift 

Now suppose that in addition to the static field, we apply a rotating field that 



is slowly varied as described in Section 4.2. We have seen that the Berry phase 
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acquired by a spin depends on its transition resonance frequency as given by 
Eq. (p6|). Therefore, at the end of a cychc evolution, the Berry phase acquired 
by the spin Sa will be different for the two possible states of spin St- Indeed, 
when spin Sb is in state 1 1), the Berry phase acquired by the spin Sa is 7+ = 
=F7r(l — COS0+), with the sign negative or positive depending on whether spin 
Sa is up or down, respectively, and 



(35) 



Similarly, when spin Sb is in state 1 1), the Berry phase acquired by the spin Sa 
is 7- = TT^i^ — cos6'_) where 



cos9- 



UJ- 



(36) 



As in the single-particle case, it is necessary to eliminate the dynamic phase 
in order to construct a purely geometric conditional phase gate. This can be 
accomplished using almost the same technique as in the single-particle case 



described in Section 4.3. In this case, however, we must apply the sequence of 
operations 



C 



C 



c 



(37) 



where tt^ and Wb are 7r-pulses applied to particl es a and 6, respectively, and C and 
C are adiabatic transformations as in Section O. If we define the differential 
Berry phase shift 



A7 = 7+ — 7- = TT 



then the net transformation, up to global phases, is given by 



v 





2iAf Q 













-UJt 



(38) 



(39) 



Thus, we have succeeded in engineering a conditional evolution since the state 
of the qubit Sb influences the phase acquired by a second qubit Sa- |2^] This 
gate, which introduces a phase of e^'^''' if the two spins are aligned and e^^"^''' 
if they are anti-aligned, is equivalent to the controlled phase gate introduced in 
Section m. 



5.3 Fault Tolerance 

The form of the dependence of A7 on the detuning LUa — iu and the amplitude 
of the oscillating field lui builds into the geometric phase gate a natural type of 
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Figure 5: Plot of differential phase shift A7 as a function of "° and 

fault tolerance not present in the simple non-geometric conditional phase gate. 
In many experiments, such as NMR, it is easy to control the detuning quite 
precisely, but relatively difficult to control the amplitude of the oscillating field. 
Figure |^ plots A7 as a function of these parameters in units of tt J for a range 
of values: for fixed — ^, we see that there is a peak in A7 as a function of Wi. 
Therefore, if uJi is chosen to coincide with this peak then A7 will be insensitive 
to errors in loi to first order. As the height of the peak depends on the detuning, 
any desired controlled phase gate can be achieved. 

6 Conclusions 

The techniques described in this paper constitute a novel approach to quan- 
tum computation, one that builds entangling gates entirely out of conditional 
geometric phases. These techniques are readily implementable with current 
technology in quantum optics and have already been demonstrated by some 
of the authors using NMR jl9| . It would be interesting to further analyse the 
robustness of geometric quantum computation to errors. While it has been ob- 
served that geometric phases are robust to certain types of noise in the classical 
parameters controlling the Hamiltonian, it has not been determined how geo- 
metric phases behave in the presence of decoherence or depolarisation of the 
quantum system. 
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